Analysis of a Tubular Gas Lens 

By D. MARCUSE and S. E. MILLER 

(Manuscript received June 8, 1964) 

// a cool gas is blown into a hot tube, it acts as a positive lens which will 
focus a light beam passing through the tube. 

Using a theory presented in Ref. 3, we give curves which show the temper- 
ature distribution in the tube as a function of the distance from the tube axis 
and also as a function of the distance along the axis. 

The focusing power of the lens is described by the difference in phase angle 
between a ray on the lube axis minus a ray at arbitrary distances from this 
axis and also as the second derivative of the phase angle on the axis of the 
tube. The phase curves, as a function of distance r from the tube axis, follow 
very closely an r 2 dependence. Expressions are given for the focal length of 
the lens. 

The power consumption of the lens is discussed, and a figure of merit is 
defined as focusing power per watt. The gas used for this lens should be se- 
lected such that (n — l)/k is as large as possible (where n is the index of 
refraction, k the heat conductivity of the gas). 

Using C0 2 and a \-inch ID tube 5 inches long heated 20°C above the in- 
coming gas, a focal length of 5 feel with a power consumption of 0.325 watt 
is calculated; the focal length is inversely proportional to power consumption 
within certain limits. 

I. INTRODUCTION 

A communications system using light as the carrier of intelligence 
needs an efficient medium to propagate light from transmitter to re- 
ceiver. Among the several alternatives, the idea of Goubau and Schwer- 
ing 1 of confining and propagating an electromagnetic wave with a system 
of lenses appears promising. However, in a lens-waveguide system there 
is a wide range of possibilities as to what types of lenses to use. Conven- 
tional glass lenses present problems, since they may not only absorb 
light in the glass medium itself, but furthermore present important re- 
flection losses which can be only partially avoided by special techniques 
such as coating the lens surfaces or making use of the Brewster angle. 
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Even if such corrective measures are used, there is still residual reflection 
and scattering of light due to unavoidable surface irregularities. 

It appears that most of the problems connected with glass lenses could 
be overcome if, instead of a high-index medium such as glass, a very low- 
index focusing medium were used. If the transition from air into a dense 
medium could be avoided, the problem of light reflection would not exist. 
Gases present themselves as an obvious choice of a low-index dielectric 
medium. Their dielectric constant can be influenced by changing their 
density. A change of density is most easily effected by varying the gas 

temperature. 

D. W. Berreman 2 built a successful gas lens by maintaining a tempera- 
ture gradient between a hot helix and a cold cylindrical enclosure. Al- 
ternatively, D. W. Berreman and S. E. Miller proposed a gas lens formed 
by blowing a cool gas into a hot tube (Fig. 1). Since the gas heats up first 
at the wall of the tube and remains cool longer at its center, it has a 
density distribution of higher-density gas in the center of the tube and 
decreasing density towards the wall. Since an increase in density is ac- 
companied by an increase in dielectric constant, it is easy to understand 
that the cool gas flowing through the hot tube acts as a positive lens and 
tends to focus a light beam traveling along the axis of the tube. 

We present in this article some theoretical results of the temperature 
distribution in the gas and the difference in phase angle between two 
light beams, one traveling along the tube axis and the other traveling 
closer to the wall of the tube. This phase difference is a measure of the 
focusing power of the lens. For an economical lens we want maximum 
phase shift with a minimum of thermal power. We present curves show- 
ing the power consumption of the lens as well as the ratio of phase differ- 
ence to power consumption. These data allow the construction of an op- 
timum lens. Different gases give different lens properties. A gas is most 
efficient if the ratio (n - \)/k is large, where n is the index of refraction 
of the gas and k is its heat conductivity. 

II. TEMPERATURE DISTRIBUTION 

The theory of temperature distribution in a cool gas which is blown into 
a hot tube of constant temperature is presented in Ref. 3. 

WARM TUBE 

COOL GAS (TEMPERATURE, T W ) WARM GAS 

IN f\ \ OUT 



z=o y 

Fig. 1 — Gas lens using forced flow in a tube. 
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It is assumed that the gas flow is laminar and has the radial velocity 
distribution of a viscous fluid 

v(r) = %U - (r/a) 2 ] (1) 

where 

r = distance from tube axis 

a = radius of tube 

v = gas velocity at r = 0. 

The temperature T of the gas is given as 

e = T w - T 

where T w is the temperature of the wall of the tube. It is normalized with 
respect to 

0o = T„ - T Q 

with To being the temperature of the cool gas before it enters the hot 
tube. 6/ d is expanded in terms of functions R n {r/a), which are shown in 
Fig. 2 for n = 0, 1 and 2. Values of R„(r/a) are listed in Table I. The 
temperature depends on the distance z measured from the beginning of 
the hot tube, the gas velocity vo , and the following material parameters 

k = heat conductivity measured in cal/cm sec deg) 

(deg = degrees Kelvin) 
p = gas density in gram/cm 
c p = specific heat at constant pressure in cal/gram. 

All these parameters depend somewhat on the temperature but are con- 
sidered constant in the derivation of the theory. They enter the equations 
in the combination 

a = k/avopCp . (2) 



Table I — R Functions of Fig. 2 



X 


/?■>(*) 


F(x) 


Ki(«) 


Rt(x) 





1 


1 


1 


1 


0.1 


0.9819 


0.9805 


0.8923 


0.753 


0.2 


0.9290 


0.9261 


0.6067 


0.206 


0.3 


0.845G 


0.8432 


0.2367 


-0.290 


0.4 


0.7382 


0.7382 


-0.1062 


-0.407 


0.5 


0.6147 


0.6175 


-0.3399 


-0.204 


0.6 


0.4833 


0.4880 


-0.4317 


0.104 


0.7 


0.3506 


0.3535 


-0.3985 


0.278 


0.8 


0.2244 


0.2244 


-0.3051 


0.278 


0.9 


0.1069 


0.1041 


-0.1637 


0.144 


1.0 
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Fig. 2 — Functions R , R\ and Rz vs r/a. 

The first three terms of the infinite series describing the temperature dis- 
tribution in the tube are 



e 



= 1.477 exp (- 7.316a ?J R (?\ - 0.810 exp f -44.36<r Z - 
■ Rl (1) + 0.385 exp(- 106. -*) ft Q =* 



(3) 
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The approximation is fairly poor at z = 0. However, the exponential 
factors in the higher terms of the series drop off very rapidly as z in- 
creases so that the approximation is already very good for values of 

<r(z/a) > 0.01. 

Fig. 3 shows 0/0o at /• = as a function of <r(z/a). It is apparent that for 
<r(z/a) > 0.05, the distribution of 0/0 o drops off exponentially. Fig. 4 
shows the r/a dependence of 0/0 o for different values of <r(z/a). 

To make Figs. 3 and 4 more meaningful, we list in Table II the material 
parameters for several gases at 20°C and a pressure of 700 mm Hg. 



III. POWER CONSUMPTION 

The principle of operation of our gas flow lens requires that we heat 
the cool gas inside the hot tube. Even if we neglect all power losses to the 
environment, we have to spend a certain amount of heat power to oper- 
ate our lens. For a subsequent study of lens efficiency we need to know 
this basic power consumption. At any given length z of the tube we ob- 
tain the power absorbed by the gas as 



'(z) = [ [T(r,z) - r ]pc p v(r)2Trrdr 
Jo 

= 2tt pc p v^ J r 1 - (-) 



e(r,z) 



(4) 



dr. 
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Fig. 3 — Normalized gas temperature on tube axis vs normalized distance 
along tube. 
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Fig. 4 — Normalized gas temperature vs radial position with longitudinal posi- 
tion as a parameter. 

In order to be able to perform the integration easily, we restrict our- 
selves to values of 

<r{z/a) > 0.05 

which allows us to express 0/0o by the first term of (3). In addition, we 
replace R (x) by 

r ( x ) &; F(x) = 1 - 2.06 x + 1.06 x\ (5) 

This approximation deviates no more than 2.5 percent from the actual 
value of R (x). The values of F(x) can be compared to those of R (x) 
in Table I. The integration can now be performed easily, and we obtain: 



P = - apCpVoda 



1 - 0.820 exp (-7.316.T -J . (6) 
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Table II — Gas Parameters vs Temperature 



Gas 


k 

(cal/cm sec 

deg) 


p 
(gram/cm 1 ) 


c p 

(cal/gram 
deg) 


avoo = k/pcp 
(cm'/sec) 


n - 1 


(« - 1)/* 

(cm sec deg/ 

ca!) 


CO 2 

NH 3 

CH 4 

Air 

H 2 

He 


3.93 10-* 
5.90 10- 6 
7.80 10" 6 
6.28 10- 5 
41.0 10~ 5 
35.0 10" 5 


1.84 10" 3 
0.72 10- J 
0.67 10- 3 
1.21 10-> 
0.084 10~ 3 
0.166 lO- 3 


0.199 

0.523 

0.528 

0.240 

3.39 

1.25 


0.107 

0.157 

0.220 

0.216 

1.44 

1.69 


4.20 10"* 
3.48 10-* 
4.13 10-* 
2.73 10-* 
1.23 10-* 
0.34 10-* 


10.7 
5.9 
5.3 
4.35 
0.30 
0.097 



Using p and c p from Table II, P is in calories/sec. Fig. 5(a) shows the 
power consumption as a function of normalized lens length, az/a. An 
alternative form of (6) brings out the dependence of P on the flow 
velocity v more clearly: 



r = \ kzd Q 7. 



*[l - 0.820 exp (-7.316 ?)]. 



The quantity 



V = kz/apCp = <r(z/a)vo 



(7) 



(8) 



has the dimension of velocity and is characteristic of the gas and the tube 
geometry. Fig. 5(b) shows the power consumption as a function of nor- 
malized gas velocity, v /V. 

The ratio of V/va can be related to the time / = z/vu which it takes the 
gas particles on the axis to traverse the tube of length 2 with the velocity 
Vq and to a time t which is defined by 



dT(0,t) 



(9) 



T„ - T(0,l) ' 



t is characteristic of the heat diffusion rate in a gas which rests in a tube 
whose wall temperature is T w . At. t = the gas has the uniform temper- 
ature To < T w . Its temperature at a given radius r and time t is T(r,t), 
so that T(Q,t) is the gas temperature at the tube axis at time t. 1/t is 
the time rate of temperature rise on the axis per degree of temperature 
difference between wall and axis. 
It is shown in Appendix A that 



V/v = az/a = 0.173(/ /t). 



(10) 



This shows that V/v expresses the ratio of the time it takes the gas par- 
ticles (on the tube axis) to flow through the tube of length 2 to the heat 
diffusion rate on the tube axis. Equation (10) may be used to replace 
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Fig. 5 — (a) Normalized power flow to the gas vs longitudinal position, (b) 
Normalized power flow to the gas vs normalized gas velocity. 

V/vo in (7) and in the other places where V/v or az/a appears, resulting 
(for example) in 

exp (-7.316 ~\ = exp (-7.316 ?) - exp (-1.265 ^ . 

The parameters lens length z and gas velocity v are important variables 
for other reasons, and the first two forms of the exponential may be pre- 
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f erred. We may note, however, 

t = 0.m(apc p /k) 

and typical values are 0.161 second and 0.08 second for C0 2 and air, re- 
spectively, when a = 0.125 inch. It is surprising that this time constant 
is so short. 

Another physical meaning one can give V is that it represents that 
velocity of gas flow along the pipe axis which assures that 0/0 o drops from 
its initial value of one at the beginning of the tube to 

0/0o = 9.910" 4 « 10 -3 

on the axis at its end. 

IV. FOCUSING ACTION 

A lens focuses because the optical path length varies for rays traveling 
at different distances from its axis. 

We describe the focusing action of our lens by the phase angle of a ray 
traveling parallel to the axis of the structure. The phase angle is given by 



<l> 



(r,z) = p [ n(r,x) dx. (11) 

J It 



Here, (So = 27r/Ao is the free-space propagation constant of the light beam, 
and n is the index of refraction of the gas. 

n(r,x) = 1 + (no - 1) =^- (12) 



T 

JO , _ "u 



T " 6 , /. o\ '" To\ V' (i:0 



*0- 



r V 0o/ 



The last step is an approximation for 0o/7 1 o « 1. The temperature in 
(12) has to be expressed in degrees Kelvin. 
We decompose $(r,z) into two parts: 

Hr,z) - p + *(r,*). (14) 

The first part 



<p = ft 



1 + (no - 1) 



(■ - a 



(15) 



is independent of the position /• of the ray in the gas lens, while the sec- 



(16) 
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ond part 

* = A(% - 1) £ f ^ * 

- ««. - ') %v h 2R ° 1 1 - exp (- 7 - 316 31 

-0.0183E 1 Q[l-exp(-44.3l)] 

+ 0.00363ft 2 (A |~1 - exp (-106 IjJ + • • • j 

accounts for the different amounts of phase shift in different parts of the 

lens. 

The difference between the phase angle of a ray traveling along the 
lens axis and the phase angle of a ray traveling at a distance r from the 
axis is 

A* = Mno - 1) ^{0.202(1 - R ) [l - exp(- 7.316 £)] 

- 0.0183(1 - fli) [l - exp (-44.3 ^jl (17) 

+ 0.00363(1 - R*) Tl - exp (-106 JjJ + ■ ■ ■ j. 

This form of A* shows clearly its dependence on flow velocity for a fixed 
tube length z. To study the dependence of A* for fixed flow rate and vary- 
ing length, the following form is preferable. 

A<J> - 0o(a/<r) («o - 1) (e /T ) {0.202(1 - Bo) 

.(! _ exp (-7.316crz/a)) - 0.0183(1 - 7?i) (1 - exp(-44.3<rz/a)) (18) 

+ 0.00363(1 - Rt) (1 - exp(-106<Tz/a)) + ••■}. 

Fig. 6 shows a plot of A<I> versus length of lens for r/a = 0.4. That 
means that we compare the phase difference between a ray on the axis 
and another at distance r = 0.4a from the axis of the lens. 

Fig. 7 shows the phase difference at a fixed length z as a function of 
gas velocity. In this case, A$ goes through a maximum which for 
r/a = 0.4 occurs at v /V = 6.9. The position of this maximum depends 
somewhat on the radius r of the ray used for phase comparison with the 
axial ray. Appendix B gives the theory and Table III gives the values of 
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Fig. 6 — Phase difference between ray at r = and ray at r = 0.4a vs longitu- 
dinal position. 



the position of the maximum vo/V for different values of r/a. Table III 
shows that the position of the maximum does not change much with the 
radial position of the reference ray. 

We can explain physically this maximum in A<l> versus Vo as follows. At 
very low gas velocity the majority of the gas in the tube is at the same 
temperature — the temperature of the walls, T w . It is heated up in a time 
r [see (9)] after entering. With little temperature difference between the 
gas at r = and at r > there is little A4>. As v increases, the gas on 
the axis remains at or near T , but that nearer the walls is heated because 
it flows more slowly [see (1)] and larger A$ develops. Beyond some ve- 
locity, further increases in velocity cause the gas at r = 0.4a (for ex- 
ample) to leave the tube at lower and lower temperatures — i.e., less 
temperature difference will exist between r = and r = 0.4a because the 
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0.14 




Fig. 7 — Phase difference between ray at r = and ray at ;• = 0.4a vs gas velocity. 

transit time of the gas through the tube becomes less than the thermal 
diffusion time constant r. Thus, at high velocities the A<& decreases. Note 
that this explanation (and the theory) depends upon laminar flow of the 
gas. If there is radial mixing of the gas, less A<£ would be expected than 
predicted above, and the maximum in A<S> versus v might not occur. 

The second derivative d$>/d(r/a) is a good measure of the effective- 
ness of the lens for light rays close to its axis. For a glass lens 

<£ = nf3od(r) 

where d(r) is the thickness of the lens as a function of the distance from 
its axis. The radius of curvature R of the glass lens is given by 

R n/3 dr* ' 

In order to be able to take the second derivative of <3? we have to express 
the functions Ro , Ri , and R 2 by power series with respect to r/a. For 



Table III — Maximum vo /V for Values of r/a 



r/a 
Vo/V 



0.2 
6.73 



0.4 
G.9 



0.6 
8.26 
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the second derivative on the axis at /• = 0, it is sufficient to know the 
coefficient of (r/a) 2 in the expansion. Jakob 3 gives a series expansion of 
the 72 -functions. 

R n = 1 - iPn (r/a) 2 ± ■■■ (19) 

with ft, = 2.705, ft = 6.66, and ft = 10.3. 
We get these values 

tfy) -^-^^M 1 — (-"■•9] 

+ 0.192 1 - exp (-106—1 + • 



The maximum of this curve as a function of Vo/V appears at Vo/V = 0.75. 

Fig. 8(a) is a plot of a normalized value of d 2 $> / d(r / a) 2 as a function 
of vo/V, while Fig. 8(b) shows it (with a different normalization) as a 
function of a(z/a). 

The r/a dependence of A<J> is shown in Fig. 9. 

In order to show what values the phase difference might actually as- 
sume, and also to compare different gases, we have plotted M> in Fig. 
10 for several gases and the following geometry and flow rate: 

2a = 0.25 inch 

Vq = 212 cm/sec, corresponding to 2 liters/minute or 4.77 miles/hr. 

ft = 1.07 X 10 ; ' cm -1 , corresponding to X = 5890A 

T w = 343°K 

T = 293°K 

and with the values of n — 1 and a as listed in Table II. These curves 
assume a tube length z so long that no further A<i> would be realized with 
a longer z (i.e., outgoing gas at uniform temperature). 

It is interesting to compare the r/a dependence of A<£ to the simple 
function c(r/a) 2 . For this purpose we use (17), which is written so that 
A<1> = at ?* = 0. Fig. 1 1 gives the normalized value of A4> as a function 
of r/a for several values of the gas velocity vq/V. For comparison the 
function c(r/a) 2 is shown by dotted lines. The constant c is adjusted so 
that both curves coincide at r/a = 0.4. The actual curves of A<£ are 
surprisingly close to the simple square law dependence in all cases. If the 
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Fig. 8 — (a) Normalized d?$/d(r/a) 2 vs gas velocity; (b) normalized d^/d{r/aY 
va longitudinal position. 



gas lens could be treated as a thin lens, it would act very similar to a 
glass lens with spherically curved surfaces. 

However, the gas lens is not thin and the question presents itself: 
how do different sections of the lens contribute to the over-all focusing 
effect? Fig. 12 shows the phase difference A* between a ray on the axis 
at r = and a ray at r for a fixed value v /V = 0.9 for different sections 
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Fig. 9 — Normalized A* vs r/a with longitudinal position as a parameter. 

of the lens. The curve showing i&> for the section —> z has already been 
shown in Fig. 11. The other curves show A* in the first £ of the lens 
[curve — > (f )z], the second ^[curve (f )z —> (§)*], and the last |[curve 
(t)« - * A- The contributions are surprisingly different at different radii 
and do not resemble simple (r/a) 2 dependences. However, they all add 
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Fig. 10 — A* vs r/a when a = 0.125 inch, ?><> = 212 cm/sec (corresponding to 
2 liters/min), X n = 5890 A, T w = 343°K, T = 293°K and aZ/a > 1. 

up to the — > 2 curve which does resemble the (r/a) 2 dependence very 
closely, as was shown in Fig. 11. 



V. FIGURE OF MERIT 

The focusing action of the gas lens becomes independent of the length 
of the lens if <r(z/a) > 1, as Figs. 6 and 8(b) show. We also know that, 
for a fixed length of the lens, there is an optimum flow velocity, as shown 
by Figs. 7 and 8. 

For practical applications one would like not only to obtain an effective 
lens but also to do so with a minimum expenditure of power. It is, there- 
fore, interesting to study the lens action, that is, d%/d{r/a) 2 , per unit of 
applied power. We may introduce the ratio 



M = i 



d% 



P \d(r/a) 2 /r=o 



(21) 
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0.005 



Fig. 11 — Normalized A-l» vs r/n with gas velocity as a parameter; dotted curve 
represents variation as (r/(i) 2 , normalized to ordinate at r/a = 0.4. 



as the figure of merit of the lens. From (6) and (20) we obtain 
2/3 Wo -if F 



M = 



■7'n k 



1 - 0.820 exp (-7.316 -) 



:22a) 



with 



F = 0.738 1 - exp (-7.316 -) 
.[l-«p(-4ul)] 



- 0.40;") 



0.192 



I - exp i -106- 



(22b) 
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Fig. 12 — A* contributed by first third, second third, and output third of gas 
lens for Vo/V = 6.9. 

For a given value of v /V, the figure of merit is proportional to 
(no — l)/k. It is advantageous to make this number as large as possible. 

The figure of merit M, given in (22b), is plotted in Fig. 13. 

The gases in Table II are arranged in decreasing order of (n — l)/fc. 
Of all the gases listed in that table, carbon dioxide is best suited for a gas 
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Fig. 13 — Figure of merit given by (22a) and focusing power (as in Fig. 8) vs 
gas velocity. 

lens. This does not mean, however, that gases with larger values of 
(no — l)//c cannot be found. 



VI. FOCAL LENGTH 



We have seen that A<I> varies nearly as (r/a) 2 . In the region where the 
lens is weak, we may treat it as a thin lens and obtain directly a simple 
expression for focal length. 

For any thin lens it may be shown that the focal length / is given by 



/ = |/3 (r 2 /A*) 



(23) 



where 

A3> is the phase shift added on axis as compared to that for a ray at 
radius r 

/3o = phase constant of the region surrounding the lens. 
We obtain A$ from (17), which is given by the following for the gas 
velocity set to maximize A<£ at (r/a) = 0.4 — i.e., at (v /V) = 6.9: 

A* = O.839(r/a)\d o /T o )0oz(n o - 1). (24) 
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Putting (24) into (23) we obtain the following expression for the focal 
length of a weak gas lens : 

/- 0.596*-' T ° n . (25) 

z 0o(n o - 1) 

If a = 0.125 inch, 2 = 5 inches, T = 293°K, and 6 = 20°C, we find 
/ = 5 feet using C0 2 as the gas and / = 8 feet using air as the gas; the 
power transferred to the gas with C0 2 would be 0.0775 cal/sec = 0.325 
watt. 

When the gas lens is not weak, one should take into account that the 
refractive index varies both with radial position and with longitudinal 
position. Work is under way to analyze this very difficult situation. A 
simpler approach, and one which should give a first-order answer for 
gas lenses operated near the velocity producing maximum A<E> [see (17)], 
is to assume a medium within the lens 



n(r,x) = ti„(1 - |a/) (26) 



where 



x = distance (within lens )from start of lens 
/• = radius 
n a = index of refraction on the axis. 

For the gas velocity (v /V) = 6.9 it may be shown that 

*-±?fi(fio-l). (27) 

n a a? To 

In other unpublished work the authors have shown that the radial posi- 
tion of a ray (or of the axis of a Gaussian beam mode) is 

- Ti — 

r = ri cos \Za 2 x + —j= sin y/a& (28) 

where 

r» = displacement of ray at lens input 
Ti = slope of ray at lens input. 

We can use this general result to specify the focal length of a strong (or 
weak) gas lens with reference to Fig. 14. All input rays with zero slope 
will converge to a point on the axis a distance d beyond the output face 
of the lens (Fig. 14), where 

d = -t= cot Wa,t), (29) 

Va, 2 
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Fig. 14 — Diagram defining focal length and position of the equivalent thin 
lens for a thick gas lens. 

/ is the length of the lens, and n a has been taken as unity. An equivalent 
thin lens may be located a distance / back from the focal point, where 



/ = 



a 2 sin {\/azt) ' 



(30) 



This expression for focal length is valid up to (\/<M) = f/2, at which 
point d = 0. For (y/a$) > x/2 the rays cross within the lens, per (28). 
For (V«2 1) « 1 it may be shown that (30) passes into (25) and the 
location of the equivalent thin lens is in the center of the distributed lens. 
E. A. J. Marcatili 4 has solved Maxwell's equations for a medium 
characterized by (2G) and has found the normal modes for a sequence of 
lenses composed of segments of such a medium. This work relates closely 
to a sequence of gas lenses described in this paper. Experiments with 
tubular thermal gas lenses are reported by A. C. Beck. 5 



VII. CONCLUSION 

When a cool gas is blown through a warmer tube, the gas at the axis 
has a lower temperature than that near the walls. Thus the density and 
refractive index is larger at the axis and a converging lens is formed. If 
the tube were at a lower temperature than the input gas, a diverging lens 
would be formed. 

There is an optimum gas velocity for maximizing the focusing power 
of such lenses, and expressions are given for this velocity. It turns out that 
the optimum transit time for gas through the tube is approximately the 
time constant for temperature changes in a gas at rest in the tube, which 
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for typical gases (air and carbon dioxide) is about 0.1 second in a J-inch 
ID tube. Although not discussed in this paper, it is found that a J-inch 
tube 6 inches long yields (at the optimum velocity for focusing power) a 
Reynolds number well below that at which turbulence is expected. 

Expressions are given for focal length and a figure of merit expressed 
as focusing power per watt of power transferred to the moving gas. 

The best gas is one with a maximum (n — l)/k, where n is the refrac- 
tive index and k is the heat conductivity. 

APPENDIX A 

Derivation of (10) 

The relation (10) 

V/v = 0.l73(t / T ) (31) 

can be derived as follows. 

The time development of a cool gas resting in a tube of wall tempera- 
ture T w can be described as follows 



./o ( W n - ) 

T - T w - 2(T W - T ) S <f x "' v "' 



[ w,Ji{w n ) 



(32) 



with 



X„ = (k/apc p )wn and Ju(w n ) = 0. 



At t = 0, (32) becomes T(r,0) = 7' , which is constant throughout the 
tube's cross section. 

As time progresses the exponents \„t become large, so that very soon 
the first term of the series is the only contributing factor. Neglecting all 
the terms except the first, we get 

dT(0,t) 
1 dt _ k 2 _- Q F (33) 

and substituting z = v Q to we get (31). Since we neglected all but the first 
term in the series, (33) represents the asymptotic value which 1/t as- 
sumes after the initial transients have died down. 

To obtain a feeling for the accuracy of the approximation involved in 
deriving ( 33 ) , we write down the ratio of the second to the first term in 
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the sum (32) forr = 0: 



W2./1O2) 



exp[-(X 2 - Xi)l) = 0.666 exp (-4.26 -J. 



This ratio is 10 2 for t/ T = 0.980 and is 10 1 for t/r = 0.45. The approxi- 
mation is excellent for times t ^ t and is quite good for t > 0.5r. 
For the special example used in Fig. 10 we get for 

CO-, : r = 0.101 sec 
air: r = O.OS^sec. 

APPENDIX B 

The Maximum of (17) 

We seek an expression for the value of v /V which brings A«p, equation 
(17), to a maximum: 

d(A#) . ( ,, do 

d(v /V) To 



• 0.202(1 - Ro) J I - exp (-7.316 -) (l + 7.316 -J 
- 0.0183(1 - Ri) <1 - exp (-44.3 -J (l + 44.3 - 



(34) 



+ 0.00363(1 - /?•>) <1 - exp (-100 -) (l + 106 -j\ . 

We set (34) equal to zero, and noting that the second and third exponen- 
tials are small, we neglect them (to be justified by the solutions thus ob- 
tained) yielding 

exp ( -7.310 - J (l + 7.310 -J 

V *A *) (35) 

= 1 _ 0.0900 £ - f\ + 0.018 , ( j ~ 1 2 \ . 

(1 — lio) (1 — /to) 

Equation (35) gives the approximate value of Vo/V at the maximum of 
A<I> for any chosen radius, r/a. Additional terms in (34) can be taken if 
more accuracy is desired, which was done in computing the v /V for 
r/a = 0.2, as given in the body of the paper. 
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